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^^ ' Global well-posedness of the time-dependent (degenerate) thermistor problem remains 

^^ , open for many years. In this paper, we solve the problem by establishing a uniform-in- 

time BMO estimate of inhomogeneous parabolic equations. Applying this estimate to 
P I , the temperature equation, we derive a BMO bound of the temperature uniform with 

.^ ' respect to time, which implies that the electric conductivity is a A2 weight. The Holder 

continuity of the electric potential is then proved by applying the De Giorgi-Nash-Moser 
estimate for degenerate elliptic equations with A2 coefficient. Uniqueness of solution is 
2 . proved based on the established regularity of the weak solution. Our results also imply 

the existence of a global classical solution when the initial and boundary data are smooth. 
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1 Introduction 

The thermistor problem refers to the heating of a conductor, with temperature-sensitive 
electric conductivity, by electric current. Let (j) be the electric potential and let E = Vcj) 
be the electric field. The electric current J is related to the electric field via J = cr(n)E, 
where cr{u) is the electric conductivity of the conductor, dependent upon the temperature 
u. The heat produced (per unit volume) by the electric current is given by Joule's law: 
E • J = cr(u)|V(/'p, and the conservation of charge is described by V • J = 0. 

Let fi denote the domain possessed by the conductor. Based on the above formulations, 
the temperature u and the electric potential (j) ^^^ governed by the equations 

— -V-{K{u)Vu)=a{u)\V4>\\ (LI) 

- V • (o-(u)V(/>) = 0, (1.2) 

for X £ Q and t > 0, where k(u) is the thermal conductivity. In this paper, we consider the 
above equations with the Dirichlet boundary/initial conditions: 

u{x,t) = g{x,t), (j){x,t) = h{x,t) for x £ dQ and t > 0, 
u{x,0) = uq{x) for X £ 0,. 

The mathematical expressions of cr{u) and k{u) depend on the materials. For some 
semiconductors, the electric resistivity p{u) = l/a{u) can be approximately expressed as [17] 

piu) = cToe^/"n, 

and the thermal conductivity k{u) can be regarded as constant (independent of it). For metal- 
lic conductors, the electric conductivity and the thermal conductivity obey the Wiedemann- 
Franz law [20j : 

where L = 2.44 x 10~^WilK~^ is the Lorentz number. In general, the electric resistivity of 
metals increases as temperature grows. At high temperatures, the electric resistivity increases 
approximately linearly with temperature: 

p{u) = po[l + a{u - ur)], 

where ur is some reference temperature and a is called the temperature coefficient of resis- 
tivity. If the temperature does not vary much, the above linear formula is often used. More 
precisely, the electric resistivity is give by the Bloch-Griineisen formula |29j : 

Pin) = piO) + A^-j j^ (,._!)(!_, -.) d.. 

where A, and n > 2 are all positive physical constants. 

For both metals and semiconductors, the electric conductivity o"(n) tends to zero as the 
temperature u grows to infinity. The elliptic equation (|1.2p is thus possibly degenerate, which 
leads to severe difficulties for the analysis of the coupled system. 



The non-degenerate assumption cji < criu) < (T2 is often used to simplify the problem. 
Mathematical analysis for such non-degenerate problem has been studied by many authors in 
the last two decades. Existence of weak solutions was studied by Antontsev and Chipot [6], 
Allegretto and Xie [1] and Cimatti [8]. With the same non-degenerate assumption, Elliott and 
Larsson [9j proved the existence of strong solutions for the 2D problem by using the energy 
method (and uniqueness follows). The 3D problem is much more difficult. To deal with the 
3D problem, one has to fully explore and make use of the coupling of the equations. The 
milestone was acheived by Yuan and Liu [251 [26] , who proved the existence of C" solutions 
for the 3D problem by using the method of Layer potentials. Yin [27] obtained the same 
result by using the techniques of Campanato spaces. Their results imply the existence of 
classical solutions when the boundary and initial data are smooth. 

Without the non-degenerate assumption, the problem becomes much more difficult. Xu 
[23] proved partial regularity of the solution, i.e. the solution is smooth in an open subset 
Z? C whose compleinent i}\D is a set of measure zero. Later Xu [23] proved existence 
of solutions with bounded temperature when the boundary potential is small enough, i.e. 
ll^llL°°(9r2x{o,T)) is small enough. Hachimi and Ammi [llj proved existence of weak solutions 
by the monotonicity-compacity method. Montesinos and Gallego |181I19] proved existence of 
"capacity solutions" by considering a new formulation with the transformation $ = a{u)V(j). 
Uniqueness of the weak solution and existence of global classical solutions remain open. 
Overall, the main difficulty of the degenerate problem is the lack of a L°° bound for the 
temperature u. 

In this paper, we overcome this difficulty by establishing a uniform-in-time BMO estimate 
for inhomogeneous parabolic equations with possibly discontinuous coefficients. Applying this 
estimate to the temperature equation, we obtain a uniform-in-time BMO bound of the tem- 
perature u, as a substitute of the L°° bound. Based on the BMO bound of the temperature, 
we further prove that the electric conductivity (t{u) is a A2 weight uniform in time. The 
Holder continuity of the electric potential (p is then proved by applying the De Giorgi-Nash- 
Moser estimate for degenerate elliptic equations with A2 coefficient. The Holder continuity of 
the temperature is proved by using the Holder continuity of the electric potential. Existence 
of a weak solution in a bounded Lipschitz domain is proved, and uniqueness of the weak 
solution is proved based on the established regularity of the solution. Our results also imply 
the existence of a global classical solution when the initial and boundary data are smooth. 

For interested readers, we refer to [3l[5][9] [T^ll5l[28] for numerical methods and numerical 
analysis of the thermistor problem. 

The rest part of this paper is organized in the following way. In Section [2] we introduce 
the notations to be used in this paper and in Section [3] we present our main results. In 
Section HJ we establish a uniform-in-time BMO estimate for the solutions of inhomogeneous 
parabolic equations, and in Section [5] we present Holder estimates of parabolic equations in 
terms of the Campanato spaces. Based on the estimates obtained in Section [4] and Section 
[H we prove global existence and uniqueness of a weak solution to the degenerate thermistor 
problem in Section [6l Conclusions are drawn in Section 7. 

2 Notations 

Before we present our main results, we define the notations to be used in this paper. 



Let n be a fixed positive integer and let Bfi{xo) denote the ball of radius R centered at the 
point xo € M". Let be a bounded Lipschitz domain in M", i.e. ri is a bounded domain in 
M" and for any y S dil, there exists a ball Buijj) such that through a rotation of coordinates 
(if necessary) , 

BR{y) r\VL = {(xi,--- ,Xn) G BR{y) : x„ > ip{xi,--- ,Xn~i)}, 

where if : M"~^ — )• M is a Lipschitz continuous function. For a bounded Lipschitz domain, 
there exists a positive constant Rq and a finite number of balls BR^{yi), Bji^{y2), •••, 
BR^{ym) such that dVt C ^T=iBR^/2{yj) a-nd through a rotation of coordinates (if necessary), 

B2Rn{yj) n = {(xi, • • • ,x„) G B2R^{yj) : x„ > ipj{xi, ■ ■ ■ ,x„„i)} 

for some Lipschitz continuous function fj : W^~^ — t- M. 

For any integer m>0, l<p<oo and < a < 1, let W"^''P{Q) and C"""'""(0) denote the 
usual Sobolev space and Holder space [l], respectively, and let C"^~^°'{il) denote the space of 
functions which belong to C"^~^°'{B) for any closed ball B C Q. Let C^~'""(0) be the subspace 
of C™"*'"(r2) consisting of functions vanishing on the boundary 9f2. 

Let \D\ denote the Lebesgue measure for any measurable subset D of M", and let Br[xq) 
denote the ball of radius R centered at the point xq G K". Let be a bounded Lipschitz 
domain in M". We say that a positive locally integrable function w defined on M" is a A2 
weight if 

\\B\Jb J\\B\Jbw{x) J 



for some positive constant C, where the supremum extends over all balls in B in M". 

For any measurable subset D of W^, we let fo = tjj\ j^, f{x)dx denote the average of / 
over D. For 1 < p < 00 and < < 1, let L^' (il) denote the Morrey space of measurable 
functions / such that 



sup f^;^ / |/(x)|Pdx) <oo, 



Bn{xQ)\^''" JBR(xo)nn 

where the supremum above extends over all balls Br{xo) with xq G O and < i? < Rq. For 
1 < p < 00 and 1 < < 00, let £^'^(0) denote the Campanato space of functions bounded 
(or vanishing for > 1) on the boundary d^, equipped with the norm 

1 



sup (^ f Ifix^dxY 

R(xo)nnV-f< ■ JBn(xo)nn J 



BR(xo)nn\'^"" JBR{xo)nn 

+ s^P [-^ I \f{x)-fBn{yo)nn\''dx'' 

where the supremum above extends over all balls with xq G d^, yo G fi and < i? < Rq, 
andweset BMO = £i'i(0). 

For any fixed T > 0, we set Of = x (0, T] and F^ = d^ x (0,T]. For any point 
(xo,to) G M""*"^, we set Q/j(xo,to) = Br{xq) x (to — -R^,io] as the parabolic cylinder centered 



at (a;o,to) of radius R. For integers tti, n > 0, < a, /3 < 1 and any open subset Q C JIt, 
let C'^~^"''^^^{Q) denote the anistropic Holder space of functions, equipped with the norm 

1,1 _ _ >r- \Dlf{x,t)-Dlf{y,s)\ 



^ \Dy{x,t)-D'lf{y,s) 

- 2^ sup 



'^' (?/,«)eQ 

and set C^iSlT) = C'^'''{TiT)- Let C™+"'"+'^(Q) denote the subspace of C'"+"'"+^(Q) with 
functions vanishing on the boundary dO,. Let C°°{Q) denote the space of functions whose 
partial derivatives up to all orders are uniformly continuous on Q. Let C"^~^"'"'^^{Qt) and 
C^i^flx) denote the space of functions which are in C"""'""'"+'^((5) and C°°((5) for any closed 
cylinder Q C ^Ti respectively. For any measurable subset Q of M""*"^ and any integrable func- 
tion / defined on Q, we let \Q\ denote the Lebesgue measure of Q and let /q = tAt /o f{^)^x 
denote the average of / over Q. Analogous to the Morrey space LP'^(O) and the Campanato 
space CP'^{Q), for 1 < p < oo we can define the parabolic Morrey space Lpara(^T) equipped 
with the norm 

and the parabolic Campanato space >Cpkra(^T) of functions vanishing on the boundary Tt, 
equipped with the norm 

+ ^^^P ( p(n+2)e / \fix)- fQRl^dx)" , 

QR{yo,so)nnT\-^^ ' JQR{yo,so)nQT J 

where the supremums above extend over all cylinders with xq G 90, yo £ ^i ^Oi^o G (0,T] 
and < i? < i?Q . 

For any Banach space X and time interval (^1,^2) C M, we denote by U'{j(t\^t'i)\X') the 
Bochner space equipped with the norm 

*2 \ p 



Lv{{t^,t2)\X) 



X 

ti 



\\f{t)\\ldt\\ l<p<oo, 



ess sup ||/(t)||x, p = 00. 

t£(tl,t2) 



The importance of the (parabolic) Morrey spaces is that LP'^{Q) translates just like 
LP'^^~ '[Q), i.e. through the transformation f{y) = f{Ry) we have 



LP'»{Br) 



^"('-')/^II/IIlM(b,)> 



just like 



LP/(T-fi)(BR) - ^ \\J IIlp/(i-'')(Bi)' 



for any ball Bji C i^. Similarly, Lpkra(i^r) translates just like L^'^^ ^\Qt)- Therefore, 
LP'^^in) and L^fra(f^T) can be used as substitute for Lp/^^-'^^Q) and Lp/(i-'')(J1t), respec- 
tively, with lower order integrability. The importance of the (parabolic) Campanato spaces 
includes: 

(1) £P'^(fi) are equivalent for ah 1 < p < oo, i.e. ^^'^fi) ^ BMO; 

(2) lil<e <{n + p)/n, then £P'^(fi) ^ C^{n) for a = n{e - l)/p. 

(3) If 1< < (n + 2+p)/{n + 2), then Cl'Linr) = Co^'^^^inT} for a = {n + 2){e- I) /p. 
These properties of the Morrey and Campanato spaces can be found in [71121]. 

In this paper, we let Cp^^p^^... ^p^ denote a generic positive constant which depends on the 
parameters Pi,P2-,- ■ ■ ,Pm- 

3 Main results 

First, we establish a uniform-in-time BMO estimate and a Holder estimate for the solution 
of the parabolic equation 

du -* 

— - V • (^V^) = V • / + /o, in nx (0,T), 

u = g on aJlx (0,r), (^'^^ 

u{x,0) = Uo{x) for X G ri, 

where $7 is a bounded Lipschitz domain in M" and A(x,t) = [Aij{x,t)]nxn is a symmetric 
positive definite measurable matrix function defined on M""*"^ such that 

n 

K-^\^f < Y, Aij{x,t)Ci^j < K\^\\ for all ^ G M" (3.2) 

holds almost everywhere for {x,t) £ W^'^^, where i^ is a positive constant. 

Theorem 3.1 (BMO and Holder estimates of parabolic equations) 

There exist positive constants C and oq € (0, 1) depending only on the elliptic constant K, 
the domain VI and the dimension n {independent ofT), such that the solution of (|3.ip satisfies 
the BMO estimate 

II^IIl°°((0T)-BMO) - C'(||/o|lLl.n/(n + 2)(Qj,) + ||/||i2,n/(n + 2)(!^j,) + 1 1 UQ 1 1 ^oo (f^) + | j^ 1 1 ^oo (p^ ) ) . 

(3.3) 
If the compatibility condition uq{x) = g{x, 0) for x € dQ is satisfied, then we have 

lkllc°''^/2(flT) - ^(ll/ollLl.("+«)/("+2)(nT) + ll/llL2,(™+2a)/(n+2)(Qy) + 1 1 ^Q 1 1 ^c (Q) + |bll C"'"/^ (Ft) ) ' 

(3.4) 
for < a < ao . 



The inequality (j3.3p is new. A similar inequality as (j3.4p was proved in [27], where 

||/o||Li,{n+<.)/(„+2)(j^^) was replaced by ||/o||i2,(„-2+2a)/(n+2)(Qy). Note that L2'('^-2+2")/("+2)(-j7y) 

translates in the same way as L^'^^~^'^''^'^~^'^'(Qt) under a scale transformation but requires 
higher integrability. 

Secondly, by applying Theorem 13.11 we prove global existence and uniqueness of a weak 
solution for the degenerate thermistor problem under the following physical hypotheses: 

(HI) The thermal conductivity is a smooth function of temperature and satisfies that 
< inf k{s) < supK(s) < oo, for any fixed r > 0. 

s>r s^j. 

(H2) The electric resistivity p{u) = l/a{u) is a smooth function of temperature such that 
for some p > there holds 

Cl,r + C2,rSP < P{s) < Cs.r + Ci^rS^ V S > r > 0, (3.5) 

where Ci^r, i = 1, • • • ,5, are some positive c onstants (possibly depending on r). 

Clearly, the hypotheses (H1)-(H2) are true for metals and some semiconductors. In par- 
ticular, the electric resistivity p{u) can be any polynomials which are positive for u > 0. The 
hypotheses (H1)-(H2) also imply that for any given r > 0, a{s) is bounded for s > r. 

Theorem 3.2 (Global well-posedness of the degenerate thermistor problem) 

Let Q be a bounded Lipschitz domain in M" (n = 2,3) and let qq > n. Assume that uq G 

VFi'9o(rj), g G L°^((0,r);VFi'9o(rj)), dtg G L'^{{0,T);Lio{n)), h G L'^{{0,T);W^'i''{n)), 

with 

min g{x,t) > 0, min uo{x) > 0, 

ix,t)erT xen 

and g{x,0) = uq{x) for x G d^. Then, under the hypothesis (H1)-(H2), the initial-boundary 
value problem (jl.ip - (jl.3p admits a unique weak solution (u, (j)) such that 



u G C°'"/2(0r) n LP((0, T); W^''i{^)), (j) G L°°((0, T); ^^^'^(O)), 

dtueLP{{o,Ty,w-^'i{n)), 

for some q>n, 0<a<l and any 1 < p < oo, in the sense that the equations 

/ / -^vdxdt+ / / K{u)Vu-Vvdxdt= / / a{u)\V(t>\^vdxdt, 
Jo Jn ^* Jo Jn Jo Jn 

/ / a{u)V(t>-V(pdxdt = 0, 
Jo Jn 

hold for any v,ip e ^^((0, T); H^ (n)). 

Note that with the regularity (j3.6p . the last equation above is equivalent to 



(3.6) 



a{u)V(p-V^pdx = 0, y^pGH^in), a.e. tG(0,r). 
Jn 



4 BMO estimate of parabolic equations 



The solution of (13. ip can be decomposed into three parts, i.e. the solution of the following 
three problems: 



3ii 

__V(^Vn) = /o, 

n = 
u(x,0) =0 



in r^ X (o,r), 

on do. X (0,r), 
for X G ri. 



(4.1) 



fill -* 

— -V-{AVu) = V-f, in Ox(0,r), 

n = 
u(a;,0) = 



on (9fl X (0,r), 
for a; G Q, 



(4.2) 



— - V • (ylVu) = 0, inJ]x(0,r), 



u = g 

u{x, 0) = uq{x) 



on (9J1 X (0,r), 
for X £ Q. 



(4.3) 



From the maximum principle and the De Giorgi-Nash-Moser estimates, we know that 
there exist positive constants C and < ckq < 1 such that the solution of (|4.3p satisfies that 

II^IIl°°{Ct) - ll5'llL°°(rT) + ll^o||L°°(f7); 

for < a < ao < 1 and T > (the second inequlaity above requires the compatability 
condition). To prove Theorem 13.11 it suffices to present estimates for the equations (j4.ip - 

The rest part of this section is organized in the following way. In Section 14. H we present 
local L^ estimates for the solution to (14.11). In Section 14^ we combine the local L^ estimates 



to derive a global BMO estimate based on the equivalence of BMO with the Campanato space 
£^'^(il). In Section [4.31 we establish the BMO estimate for (j4.2p in terms of the Campanato 
space £^'^($7). 

4.1 Local L^ estimates 

In this subsection, we present local L^ estimates for the solution of (j4.ip . The estimates 
obtained in this subsection will be used in Section 14.21 to derive a global BMO estimate 
uniformly with respect to time. 

Lemma 4.1 Let xq £ Q and < to < T. There exists oq E (0, 1) and C > such that if u 
is the solution of (14. ip in Q^ = B^{xq) x //j with Ir = (to — B^^t^], then 



max|K-.Qj|,.(^^)<c(|) 



n+ao 



max \\u 
taiR 



holds for any < p < R < min(dist(3;o, 5(7), -v/to) CLi^'d any 
aQ depend only on K and n. 



L^{Bn)+C\\k\\L\QR) 

G M, where the constants C and 



Proof First, we prove the lemma for 9 = 0. Let B,. = Br{0), Ir = {—r'^, 0] and F^ = dB^ x I^- 
With any function ^ defined on Qr, we associate a function CiVis) = ^(xq + Ry,to + R'^s) 
defined on Qi := Bi x Ii. Then n is a solution to the equation 

^ - V, • (AVyu) = R'fo 

in Qi. Let w be the solution of 

^-Vy {AVyw) = i?7o 

with the boundary/initial condition w = on the parabolic boundary dpQi and let w be the 
solution of 

^-Vy{AVyw) = R'mQ^ 
in M""''^ with the initial condition w{y,0) = 0. By the maximum principle, we know that 

\w{y,s)\ < \w{y,s)\ 

Taking the L^{Bi) norm with respect to y, we derive that 

II^IIl°°(7i;L1(5i)) - ^^ II-/'oIIl1(Qi)- 

We note that v = u — u;^ — w is the solution of 



dv 
'd's 



Vy ■ {AS/yV) = 



in Qi, and by the De Giorgi-Nash estimates of parabolic equations we know that there exists 
ao G (0, 1) such that for p G (0, 1/2], 



- — / \v — vp^ \dy 



max 



Therefore, 



te/i 



max ki - u^ Wtkr \ 

< m|x 11^ - ^;QJIil(B^) + m|x \\w - WqJ\^,^^^^ 

< Cp""*""" max \\v\\ r-itS \ + Cmax \\w\\ r-i/s \ 

- F ^^^^ IIL (Bi) ^^j^ ML (Bi) 

< Cp""''"'' max llu — Up; II ri,B ^ + C max ||w|| f 1,5 X 



< Cp"+"o max \\u - u^^ H^.^^^^ + CR^WM^.^q 



teh 



LHQi) 



< Cp-+-- max 11^11^,(5^) + CR^foW^.^Q^y 
tall 



where we have noted that 

„~ „ l^il f ,~, , 1 

Transforming back to the (rr, t)-coordinates, we complete the proof of the Lemma, for = Q. 
Then we note that u — \s also a solution to the equation (j4.ip in Qn for any ^ S K. D 
Similarly, we can prove the following local L^ estimates near the boundary dpVtj'- 

Lemma 4.2 Let xq G ^ and Iq = 0. There exists oq £ (0, 1) and C > such that if u is 
the solution of (j4.ip in Qr = Br{xq) x /^ with /^ = [0, i?^], then 

max||u||ii(B^) < C\^-j max||n||ii(B^) +C||/o||li(q^) 

holds for any < p < i? < min(dist(2;o, dQ), vT), where the constants C and oq depend only 
on K and n. 

Lemma 4.3 Let xq G dVl and to > 0. There exists ao G (0, 1) and C > such that if u is 
the solution of (j4.ip in Qr = Br x Ir, with Br = Br{xo) Ci Q, and Ir = (to — R^,tQ\, then 

/ pY+'^o 
max||w|Ui(B^) <C7(^-j max ||n||ii(B^) + C7||/o|Ui(q^) 

holds for any < p < R < min(i?Q, y/to), where the constants C and ao depend only on K , 
n and il. 

Lemma 4.4 Let xq G dVl and to = 0. There exists Uq G (0, 1) and C > such that if u is 
the solution of (j4.ip in Qr = Br x Ir, with Br = Br{xq) n $7 and Ir = [0, R'^], then 

f p\ "+"" 
max||n||ii(B^) < Cy^-j max ||n||ii(B^) + C7||/o||ii(Q^), 

holds for any Q < p < R < min(i?f7, vT), where the constants C and ao depend only on K, 
n and i7. 

The following simple lemma can be found in [TllTB], which is widely used for estimates in 
terms of the Morrey and Campanato spaces. 

Lemma 4.5 Let (/?(•) he a nonnegative and nondecreasing function defined on (0, i?o] and 
suppose that for any < p < R < Ro, 

where Ci , 71 and 72 are nonnegative constants such that < 72 < 71 . Then 

10 



From the above lemmas, we obtain the fohowing local L^ estimates. 

Proposition 4.6 For xq G il, to > and Qr = Bji(xq) x Iji with Ir = {to — i?^,fo]. we 
have 

for any < p < R < min (dist(xo, 50), -v/to) . 

Proposition 4.7 For xq ^ Q, to = and Qr = Br{xq) x [0, R^], we have 

-;rlkllL-((fo-p2,to);ii{Sp)) - C'( ;^ll^llL-{(to-«^to);ii(Bfl)) + ll/ollLi."/("+2){nT) 

for any < p < R < min (dist(xo, dU), vT) . 

Proposition 4.8 For xq G 50, to > and Qjj = Br{xo) HQ, x Ir with Ir = (to — R^^t^], 
we have 

for any < p < R < min(i?n, -v/to). 

Proposition 4.9 For xq G dQ, to = and Qr = Br{xq) n O x [0, R'^], we have 

1 „ „ 



D,p2);Li(Bp)) < C{—\\u\\La.^^Q^R2y^L\Ba)) + II/o|Ili."/("+2)(Qj,) j 

/or any < p < R < min(i?Q, vT). 



^ll«llL-((0,r2vri^ 



4.2 BMO estimates via £i^ 

We combine the local L^ estimates obtained in the last subsection to derive a global BMO 
estimate of u, uniform with respect to time. 

Proposition 4.10 The Propositions I4.6fl4.9l imply that the solution of (|4.ip satisfies that 

II^IIl°°((0,T);BMO) - <^ll/ollLi^"/("+2)(nT)' ('^•'^) 

where C depends only on K, n and O (independent ofT). 

Proof Set M = ||/o||ii,n/(n+2)(!^j,). 

First, we prove the proposition for T > Rq. We shall prove that for R < Rn/2 and any 
set Br = Br(xq) n with some point xq G and 5 = dist(a;o, 90), the following estimates 
hold: 

J w\\'"-\\L°°i(0,Ty,L^iBit)) - ^{\W\\L°°{iO,T);L^{n))+ M), if 6 < R, 

[ 7rr||u - ■UB^||ioo((o,T);Li(Bfl)) - ^(ll'"lli-°°((0,T);Li(n)) +^)) if S > R, 



11 



Case 1: 6 < R. In this case, there exists a region B2R = B2R{yo) H $7 with some yo S 50 
such that Br C i32R and so, for any given to ^ i^jT], 



l^(->io)||Li(B«) < ll'"(-,io)||Li(B2fl)- 



(4.6) 



Now if io < 4-R^, then by Proposition W^ 



;^lh(->*0)||Li(B2fl) < lhllL°°((0,4ij2);Ll(B2fl)) 



1 



< C lln 



j^n\m\L-°{iO,Rly,LHBR^^)) 



+ M 



(4.7) 



Otherwise, to > 4i?^ and by Proposition l4.8l for Rq = min(-y/to, Rq) and Rm = max(T/to, Rq) 
we have 

1 






< < 



C( 7ir||n||ioo((o,R„);ii(B^^)) +M ), 



/to < ^n, 

/t^ < Rn, (by Proposition liT9 

C( 7|^ll^llL°°({to-fl^,io);iHBfln)) + ^'^ ) ' ^-^^■ 
To conclude, for 6 < R and to G [0, T] we have 
1 C 

-^\\u{-,to)\\mBn) < ^^ll'"llL-({0,T);Ll(n)) +<^Af. 



i? 



Rr^ 



{4.i 



Case 2: 5 > R. In this case, we set Rq = m.m(6,y/to,Rfi). Then Proposition 14.61 imphes 
that 



1 



1 



-^,\\U - UBii\\L--i{to-R2,to);L^Bn)) < ^\\U - MQfl||L-((io-R2,to);Li(Bfl)) 

C( 7|fll"llL-((io-R?i,io);ii(Bfln)) +^^)' ^^ ^^ - min(<5,\/to) 
C( ^l|w||L-((to-52,to);iHS5)) +^]' else if 5 < min(^, i?n) 



C| 7|"ll"llL°°((0,i?2);Li(Bflo)) + ^^ 



else to = Rq 
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wII«IIl-{(0,T);L1{O))+C'^, 



< < 



c_ 

B^\W\\l°-({0,T);L^(Q)) 



+ CM, 



C{ ^||n||ioo((o^52).^l(5_j)) + M 



if Rn < mill ((5, y/t^) 

else if 6 < min(^, Rq) by (gS 

else if 6 < Rfi 



Ci 7^||'"||i^((o,ij2^);Li(B^^)) + M ) , else if 6 > Rq 



(by Proposition 14.91 



C7 



^ ^^II^IU°°((o,T);Li(n)) + C'Af, again by (gSl). 

So far we have proved (j4.5p . Once we note that ||^i||L°°((o.T);Li(r2)) ^ C'||/o||Li(nj,), we 
derive (fO|) from ([i3]) . 

Secondly, we prove the proposition for < T < i^f^. In this case, we consider the solution 
u of the equation 



— - V • (AVu) = /o 



(4.9) 



in the domain fi/jj^ = Q x (0,Rq) with the boundary and initial conditions u = on dQ x 
(0, Rn) and 'u(x, 0) = for x G il, where 



fo{x,t) 



foix,t), for tG(0,r), 
0, for tG {T,Rn). 



Check that 



||/0|lLl."/("+2){nflj^) ^ C'||/o||j^l,n/(n + 2)(Q^), 

ll/o||Li(f^iJn) - C'll/ollLi(f^T)' 
II^IIl°°{(0,T);BMO) - ll^llL°°((0,i?n);BMO)' 

where the constant C does not depend on T (as T — )• 0) . Then we apply the inequality (14. 4p 
to u with T = i?f7. n 

4.3 BMO estimates via £2,^ 

In this section, we present estimates for the solution of (j4.2p . The idea is similar as Section 
14.21 From the proof of the following lemma we can see the main difference between the 
current subsection and the last subsection. 

Lemma 4.8 Let xq G fi and < to < T. There exists uq G (0, 1) and C > such that if u 
is the solution to (|4.2p in Qr = Br^xq) x Ir with Ir = (to — R'^,to], then 

max \\u - nQj|2,(^^^ < C{-\ max ||n - e\\l,^j,^^ + C||/|li2(Q^) 



te/fl 



holds for any < p < R < min(dist(xo,5r2), -^/to) and any G M, where C depends only on 
K and n. 
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Proof Let Br = -Br(O), Ir = (— r^jO] and T^ = dBr x ly With any function w defined on 
Qr, we associate a function ^{y, s) = ^{xq + i?y, to + R^s) defined on Qi := Bi x Ii. Then tt 
is a solution of the equation 



Vy ■ {AVyU) = RVy ' f 



du 
ds 

in Qi. Let w be the solution of 

^ - V, • (AVyW) = RVy ■ f 

with the initial and boundary condition w = on the parabolic boundary dpQi. Multiplying 
the above equation by w and integrating the result over Qi, we obtain that 

On the other hand, we observe that v = u — up; — tf is the solution of 

dv 

— -Vy{AVyV)=0 

in Qi. By the De Giorgi~Nash estimates of parabolic equations, we know that there exists 
ao G (0, 1) such that for p G (0, 1/2], 

Therefore, 

max||i2-'Lt^ |p ~ 

< Cmaxlb - t;;^ |P,,,5 , + Cmax llu; - t(;^ llr2/5 n 

te/p Qpi'^-^CSp) te7p Qp"i2(Sp) 

< C/9"+2°o max |b||2, 5 , + Cmax||w||2 ,5 , 

^^^"^'""-f 11-11^2(5.) + ^^'ll^~lli^(Qi)- 

Transforming back to the {x, t)-coordinates, we complete the proof of the Lemma for 9 = 0. 
Then we note that n — is also a solution to the equation (|4.2p in Qn for any G R. D 

In a similar way, we can prove the following lemmas and propositions. 
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Lemma 4.9 Let xq € ft and to = 0. There exists oq S (0, 1) and C > such that if u is 
the solution of (14. 2p in Qr = Br{xq) x /^ with /^j = [0, i?^], then 



n+2ao 



max ||u||i2(5^) < C(^-J max HnUi^^^^^ + C\\f\\l,^^Q^^ 
holds for any < p < R < min(dist(3;o, d^), vT), where C and ao depend only on K and 



n. 



Lemma 4.10 Let xq E 90 and to > 0. There exists ao E (0, 1) and C > such that if u is 
the solution of (j4.2p in Qr = Br x Ir with Br = Br{xo) n 0, and Ir = (to — -R^, to]? then 



max||n||i2(s^) <C(^-J um^\\u\\l,^j,^^+ C\\f\\l,^Q^^ 

holds for any < p < R < min(i?n, y/to), where C and ao depend only on K, n and O. 

Lemma 4.11 Let xq E di} and to = 0. There exists oq E (0, 1) and C > such that if u is 
the solution to (j4.2p in Qr = Br x I_r, with Br = Br{xq) n Q and I_r = [0,-R^], then 

m|^ \MhiB,) < C(^;^J max IMl^f^R^^ + C\\f Wl^^^^^ 



holds for any < p < R < min(i?Q, VT), where C and ao depend only on K , n and Q,. 

From the above lemmas, using Lemma 14.51 we can derive the foUowing results concerning 
the solution of (lOD. 



Proposition 4.12 For xq E fi, to > 0, Qr = Br{xq) x (to - -R^,to] and < p < R < 
min (dist(xo,5ri), "v/to), we have 

h - '^Qp\\L'^{{to-p^,toy,LHBp)) ^ ^( ;^ll"llL°°((io-H2,to);L2(i3fl)) + ll/llL2,n/(„+2)(j^j,) J P*^- 

Proposition 4.13 For xq € Q, to = 0, Qr = Br{xo) x [0, i?^] and < p < R < 
min (dist(xo, dQ), Vt), we have 

Proposition 4.14 For xq E dQ, to > 0, Qr = Br{xo) n $7 x (to - -R^, to] and < p < R< 
m.m{RQ,^/to), we have 

ll"llL°°{(to-p2,to);L2(Bp)) < C"! ;^ll^llL°°((io-R2,fo);L2(Bfl)) + ll/llL2,n/(n+2)(nT) ) ^"• 

Proposition 4.15 For xq E dQ, to = 0, Qr = Br{xq) r\9. x [0,i?^] and < p < R < 

min(i?f7, vT), we have 



^llL°°((0,p2);L2(Bp)) ^ C"! ;^ll^llL-'((0,ij2);L2(Bfl)) + ll/llL2,n/(n+2)(t^j,) ) P 
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With the above propositions and following the outline of Section 14.21 we can prove the 
global BMO estimate below. 

Proposition 4.16 The Propogi^iong 14. 12VI4. 151 imply that the solution of (14.21) satisfies that 

II"IIl'^((0,T);BMO) - C'||/||L2,n/(n + 2)(Qj,), (4.1) 

where C depends only on K, n and VL {independent ofT). 

5 Holder estimate of parabolic equations 

In this section, we list the propositions to be used in deriving (13. 4p . We omit the proof 
of these propositions, since it is very similar as the last section. The reason we keep these 
propositions in this section is that some of them are also used in the next section to prove 
global well-posedness of the degenerate thermistor problem. 

There exist positive constants uq and C such that the following propositions hold. 

Proposition 5.1 For xq e Q., tQ > 0, Qr = Br{xq) x (to - -R^,to], < 2p < R < 
min (dist(xo, dO,), \/to) o.i^'d < q < ao, the solution of (j4.ip satisfies that 



pn+2+ 



laT^lk - ^QpIIli(Qp) < C'(^^„+2+all^ - ^WlHQr) + -^^\\fo\\LHQn)j > 



where 9 is an arbitrary constant. 

Proposition 5.2 Forxo G ^, to = 0,Qr = Br{xo)x[0, R'^], < p< R< min (dist(xo, 90), VT) 
and < a < ao, the solution of (|4.ip satisfies that 



^n+2+a\m\LHQp) 



- ^yjln+2+aMLHQR) + -^^IT^W Ml^Q r) j ■ 



Proposition 5.3 For xq G dQ, to > 0, Qr = Br{xq) n f] x (to - -R^,to]; < p < i? < 
min(i?n, i/to) and < a < ao? ihe solution of (j4.ip satisfies that 

^n+2+all^lliMQp) - ^y^n+2+J\'^\\LKQR) + ;^^T^ ll/o||Li(Qfl) J • 

Proposition 5.4 For xq E (90, to = Q, Qr = Br{xo) n O x [0,R'^], < p < R < 
min(i?f7, vT) and < a < oq, the solution of (14. ip satisfies that 

pn+2+J\^\\LHQp) ^ ^[^Jin+2+J\^\\LHQR) + -^^^[^\\fo\\L^Qu)J- 

With the above propositions and following the outline of Section 14.21 we can derive the 
following estimate in terms of the Campanato space. 

Proposition 5.5 The solution of (14. Ih satisfies that 
where C depends only on K, n and O {independent ofT). 
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9 

The local and global estimates in £pkra(^T) follow in a similar way. To conclude, we have 

Proposition 5.6 For xq ^ Q., Iq > 0, Qr = Bn{xo) x (to - R'^,to], < 2p < R < 
min (dist(xo, dQ), -v/to) and < a < oq, the solution of (\A.2h satisfies that 

where is an arbitrary constant. 

Proposition 5.7 For xq G 0, to = 0, Q^ = -BR(a;o) n x [0,i?2]^ < p < i? < 
min (dist(a;o, 90), -v/T) and < a < ao; t^^ solution of (|4.2p satisfies that 



P" 



1 „ „2 / 1 „ „2 1 



T2T2^ll"llL2(Qp) ^ ^l ^„+2+2all"lli2(Qfl) "^ Jln+2a W J W L^Q r) 



Proposition 5.8 For xq G dQ, to > 0, Qr = Br{xo) n 17 x (to - i?^,to], < p < R < 
uim{RQ, \/to) and < a < ao, the solution of (|4.2p satisfies that 

pn+2+2all"lliHQp) - '^ ( ^n+2+2a 11^11^2(0^) + j^„+2a II/IIl2(Qb) ) ■ 

Proposition 5.9 For xo G <917, to = 0, Qr = Br{xq) r\ Q. x [0,R^], Q < p < R < 
min(i?f7, vT) and < a < ao, t/ie solution of (I4.2p satisfies that 



Proposition 5.10 The solution of ()4.2p satisfies that 

where C depends only on K, n and Vt {independent ofT). 

Proposition 15.51 and Proposition 15.101 imply the global Holder estimate p.4p . 

6 The degenerate thermistor problem 

In this section, we prove Theorem 13.21 concerning global well-posedness of the degenerate 
thermistor problem. Before we prove the theorem, we introduce some lemmas to be used. 

6.1 Preliminaries 

Lemma 6.1 Let p>0. If u e BMO{W), it > 0, and d + C2\s\p < p{s) < C3 + C4\s\p for 
s > 0, then p{u) is a A^ weight in the sense that 



for any ball B C M", where the constant C depends on Ci, C2, C3, C4,p and ||n||BMO- 
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Proof For any ball B C ffi", we set Bi = {x e B\ \u{x) - ub\ < ^ub} and B2 = B\Bi. By 
the Nirenberg inequality [H] we have I-B2I/I-BI < e-^"s/li«liBMo. Clearly, p(n) > Cp{ub) on 
B\. Therefore, 

—7 / p{u)dx<-— {l + \u-UB\^)dx + -— / \uB\^dx 
\B\ Jb \B\ Jb \B\ Jb 

<C + C\ub\p<Cp{ub), 

1 /■ 1 , 1 /■ 1 , C\B2\ 

' -dx<-— — -dx+ ' ' 



B\ Jb p{u) - \B\ Jb, p{u) \B\ 

c c 

< _j_ g-CMs/||M||BMO < 



p{ub) p{ub)' 

The last two inequalities imply that p{u) is a A2 weight. D 

The following lemma concerns maximal regularity of parabolic equations, which is an 
application of the maximal regularity of [22] and [13] (with the perturbation method for the 
treatment of operators with merely continuous coefficients). 

Lemma 6.2 Let u be the solution of the parabolic problem (13. ip in M" (n = 2,3) with the 
Dirichlet boundary /initial conditions u = g = fo = 0, and assume that the coefficient matrix 
A is continuous. Then we have 

\\u\\LP{I;W^'i{n)) < Cp^q\\f\\LP{I;Li{n)) 

for some q > n and any 1 < p < 00. The constant Cp^q depends only on p,q, K , the domain 
Q and the modulo of continuity of A. 

The analogus result for elliptic equations is given below, which can be proved by applying 
the W^''^ estimate of [13] with a perturbation argument. 

Lemma 6.3 Let Aij, i,j = 1, • • • ,n, be continuous functions defined on Q, satisfying 

n 

K-^\^\^ < Y^ Aij{x)^i^j < K\^\^, for all ^ G M", a.e. x € M" (n = 2,3), 

where K is a positive constant. Let u be the solution of the elliptic equation 

-V • (AVu) = V • / in n, 
with the Dirichlet boundary /initial conditions u = on di}. Then we have 

\W\\w^'i{n) < C'qll/llL9(n) 

for some q > n. The constant Cq depends only on q, A, the domain Vt and the modulo of 
continuity of A. 

The following lemma is concerned with Holder estimates for inhomogeneous parabolic 
equations [2] , which is also a consequence of Theorem 13.11 



18 



Lemma 6.4 The solution of (j3.ip with uq = g = satisfies that 

\\u\\(.a,c/2(Q^) < C{\\fo\\LP{{0,T);Li{n)) + ll/llL2p((o,T);L29(n))), 

for some < a < 1, provided 1 < p, (/ < oo and 2/p + n/q < 2. 

The following lemma concerns an estimate of Vu in the Morrey space for the parabolic 
equation (|3.ip . which was proved in [27] for uq = g = fo = 0. 



Lemma 6.5 The solution of ()3.ip with /q = satisfies that 



|Vn|| 2,n/{n+2).„ . < C( 1 1 / 1 1 2,n/(n + 2) ,„ s + 1 1 V^ 1 1 2,n/{n+2) , 



6.2 Construction of approximating solutions 



+ 



o?£/'"+''(f^T) + ll''oll^°°(^))- 



For the non-degenerate problem, the existence of a C" solution was proved by Yuan and Lin 
[251 I26j . Based on their result, for any given e > 0, there exists a weak solution {u^^cj)^) 
such that (/)" G L°°((0,r); //^(Jl)) and -u" G C"'"/2(nr) n L2((0,r); //^(rj)), to the following 
equations 



-^ - V • (K(n^)Vu^) = V • [(a(n") + e)4>^V4>^] in 0, 



u" = g 



u'^{x,0) = uq{x) 



on (9r2, 

for X G ri. 



(6.1) 



J -V • ((cj(n=) + e)V(/)^) =0 in 0, 
1 (^^ = /i on (9f]. 

We also note that, by the maximum principle, the solution u^ of (I6.ip satisfies that 

u^ > c := min ( minno(a;), min g{x)) > 0, 

x£f2 x&dQ 



and the solution (/) of (|6.2p satisfies that 

ll'/'^llL°°(nT) < ll^llL°°(rT)- 
By the hypotheses (H1)-(H2), we have 

kq < k{u^) < Ki, £ < (j{u^) + e < 2(To := sup a{s) 

s>c 



(6.2) 

(6.3) 

(6.4) 
(6.5) 



for some positive constants kq, ki and (Tq, where we choose e < (Tq. 
Proposition 6.1 The solution {u^,(j)^) of (j6.ip - (|6.2p satisfies that 

||w'^|lca,c</2(Qj,) + \\u'^\\LP{{0,Ty,W^-i(n)) + ll<9iM^||LP((o,r);H/-i''?(n)) + W\\L^i{0,T);W^-iin)) < C") 

and 

ll<^^llc°='"/2(BflX[0,T]) - '^dist(BH,9n) 

/or any closed ball B^ C Jl, where the constants C and CdistfU dn) ^^^ independent of e. 
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Proof First, we show that a{u^) + e is a ^42 weight, uniformly with respect to time and e. 
Let xq € ri, to > and let Rq = ^ min(-v/to, dist(xo, Sfi)). For any ball Bji of radius R 
centered at xq, we let (" be a smooth function defined on M."" which satisfies 0<C<1,(" = 1 
in Bji and C = outside i?2_R- For any interval Ir = (to — -R^jto], we let x be a smooth 
function defined on M which satisfies < x ^ 1) X = 1 oii -^-R &iid x = on (— oo, tg — 4i?^]. 
Let Qr = BrX Ir so that (M^ (/>'') is a solution of ([61])- ([621) in Q2Bo- Multiplying ([62]) by 
(f = (j)'^ C'^ , we obtain 



Jbr Jb2r 

Integrating the above inequality with respect to time and using (|6.4p . we get 

[[ ia{u') + e)\Vcl)'\^dxdt<C\\h\\l^,^.R^. (6.6) 

JjQn 

Similarly, for xq G dU, to > 0, R < ^ min(^/t^, Rq), Br := BR{xo)riO, and Qr = (5i?(xo, to)n 
Qt, we also have (16. 6p . From the last inequality we see that 



||Va(w-)+eV</.^|| 2.„/(„+2) < C. (6.7) 

By Theorem 13. H the solution of (j6.ip satisfies that 

ll^^llL°°((o,r);BMO) - C'||v/<7(m^) + eV0^||^2^n/(n+2)^^^^ + C'||^o||l°°(c) +C'll5'llL°°(rT) ^ C*- 

(6.8) 

Applying Lemma [631 to the equation ()6.ip and using ()6.7p . we derive that 

||V?i^|L2,„/(n+2),„ . < C||\/cj(u^) + eV<^"|L2,n/("+2).o , + C <C. (6.9) 

-^para V^'T; " " -t^para {'■^T) 

We extend the function u'^ defined on to M" by setting u^{x) = c for x G ]R"\ri so that 

II^^IIl°°((0,T);BMO(IR")) ^ C. 

Since (|3.5p holds, from Lemma |6. II we see that p(u^) (and also a{u'^) = l/p(n^)) is a ^2 
weight uniform with respect to time and e. It follows that, for any ball B C M", 



— - / (aiu") + e)dx] ( —- / — — ^ dx 



" VI^Ub ^ ' )\\B\jRa{u-) ) 
which says that cr{u^) + e is also a A2 weight, uniform with respect to time and e. 



20 



Secondly, we estimate the Holder norms of (j)^ and u^ , respectively. In fact, from [10] we 
know that any solution of the elliptic equation ()6.2p with the A2 coefficient (j{u^) +e satisfies 
the Holder estimates: 



{■Mc<^m<CW{-,t)\\c'^{dn)<C, fortG(0,T), VaG(0,ao), 



(6.10) 



for some fixed constant a^ G (0, 1). 

We proceed to the Holder estimate of u'^ . For any fixed xq € ri, we decompose the function 
u^ as n^ = uf + ul , where u\ and ul are weak solutions of the equations 

^ - V • (K(n^)Vn^) = in f]. 



u\= g 



nf(x,0) = Uo(x) 



on (9il, 

for X G ri, 



and 



r du\ 



-^ - V • {k{u')VuI) = V • [(0^ - 0^(xo, t))(a(n^) + e)V0^] in n, 



dt 
_ n|(x,0)=0 
respectively. By the De Giorgi-Nash-Moser estimates, we have 



on 5il, 

for X G fi, 



Ui 



lc«."/2(Qr) - C'(||5llc"."/2(rT) "*" ll'"'o|lc«(n))- 



and in order to estimate 11^*1 llca-^/^fo V ^^ ^^^ / 



(0^ - (/)^(xo,t))(cr(n^) + e)V0'' and 
apply Proposition l5.6l -Proposition I5.9i We see that for xq G J^, to > 0, < 2p < i? < 
min (dist(xo, dO), y/to), we have 



1 



,n+2+2a 



<c 
<c 
<c 



m - (^2)QpIIl2(Qp) 



1 



+ 



1 



^n+2+2Qll"2 '^llL2(Q^) -r j^n+2a\\J\\L^{QR) )^ 
1 „ . „„o 1 



^n+2+2o 

1 

^n+2+2o 



|n|-0| 



+ 



/,ei|2 



i^(Qi?,) "^ i^ri"^ llL°°(/;C"(n)) 



^rmJ|vM^^^) + eV0^ 



li2(Qfl) 



F2 



i^^Wh) ^ ;^ 



+ ]^llVcx(n^)+^v</.i^,(Q^)), 



Similarly, for xq G f^, to = 0, Qij = Br{xq) x [0, i?^] and < p < -R < min (dist(xo, Sfi), a/t) , 
we have 



1 



,n+2+2a I 



P' 



1 



1 



Hli2(Q^) < C[-—^^\\ul\\l,^Q^^ + — ||V<T(n^) + ^V0' 



i?" 



li^CQfl) /■ 



For Xq G 317, to > 0, Qij = Br{xq) r\Q.x {Iq- R^.Iq] and < p < i? < dist(i?n, \/to), we 
have 



1 



,£||2 



P' 



,„+2+2all"2llL2(Q^) 



< c 



1 



,e||2 



+ 



1 



^n+2+2all"2llL2(Qfl) "^ _^„ 
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|V<^(u^) + eV(/)^ 



\l^{Qr) ■ 



For xo G d^, to = 0, Qr = -Br(xo) n J] x [0, R'^] and < p < i? < mm{Rn, y/T), we have 



1 ....,,2 - ^/ 1 



pn+2+2all^2llL2(Q^) 



^^(;R^^K2lli.(Q,) + ;^||V^("^) + ^V0i^,(Q^)). 



Combining the last four inequahties and fohowing the outhne of Section STJl we can derive 
that 



||w||L2,l + 2a/(n + 2),„ N < C* U/(T(n^ ) + eV<?i>^ 2,n/(n + 2) , „ .. 

*-'para K^^T J " " -t^para \WR) 

With ([EZD and the equivalence relation £part^"^^"^^^(^r) = ^"'"/^(Q^^), we see that 

ll'"2llc"'"/2{nT) - ^• 

Therefore, 

lk^llc"'"/2(nT) - II'"'iIIc"'"/2(IIt) + ll'"'2llc"'"/2{nT) - ^- (6-11) 

Thirdly, we present W^''^ estimates of (p^ and u^. Note that the last inequality implies 
that 

C7-1 < a{u') + e<C, \\a{u') + e|bc.,./2(n^) < C, \\k{u')\\c...^,^j^^^ < C. (6.12) 

With the Holder estimates of a{u^) + e and ^(n^), we apply Lemma [6T2] - 16.31 and derive that 

II<^^IIl°°((0,T);W'1.9(Q)) < C'||/l||^oo((o,T);Vyi>'J(Q)) < C, (6.13) 

ll^^llLP((0,T);Vyi.'3(n)) < C'p||<A^||LP((0,r);W'i^9(Q)) + Cp < Cp, (6-14) 

for some q > n and any 1 < p < oo. From the equation ()6.ip we also see that 

l|c^t^^llLP(7;Ty-i>'!(f7)) < C'(ll^^llLP((0,T);Tyi>'J(n)) + II ^0^ ||ip((o,r);VFl''J(n))) < C. (6.15) 

Finally, we estimate the interior space-time Holder norm of (j)'^, which is used to ob- 
tain pointwise convergence of the approximating solutions in the next subsection. For the 
simplicity of notations, we set A^ = cr{u^) + e. From (j6.2p we see that 

-V • (a%x, ti)V[(/.^(x, ti) - <p%x, t2)]) = V • (iA%x, h) - A%x, t2))Vcj)%x, t^) 



By applying the interior W^''^ estimate to the above equation, we find that for any closed 
ball Bfi contained in O there holds 

\\(t)^{x,ti) - <t>^{xM)\\L'^{{0,T);W^-i{BR)) < Cdist(Bfl,af7) P^(^' *l) " ^^l^^, ^2) IIl-CCt) 

■i' n — II zi^ii _ 1+ + i"/2 

- '-'dist(Bii,en)ll^ llc"'"/2(nT)l^i ''Sl ' 

which reduces to 

\W\\c°'/'^{[0,T];W^-i{Br)) ^ ^dist(Bii,dQ.)- 

Since W^''^{Br) ^-)- C°(n), the last inequality implie that 

||0^|lc<^,a/2(;B^x[O,T]) - ^dist(BH,9n)- (6.16) 

The proof of Proposition 16.11 is complete. D 
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6.3 Existence of solution 

Since C"'°/^(r2T) is compactly embedded into C{i}T) and C"'"/^(-BrX [0, T]) is compactly em- 
bedded into C(Brx[0, T]) there exist functions u E C°''°'/'^{nT), </> S L~(/; W'^'i(n)) with (/> e 
(7":"/2(^^^ X [0,T]) for any closed ball Bji contained in Q, and a sequence e^ — ;• 0, such that 
u^'' converges to u in the norm of C{^t), u^*" converges weakly to u in LP{I; W^''^{Q)), dtu^'' 
converges weakly to dtu in LP{r,W^^''^{Q)), (f)'^'' converges weakly* to (f) in L°^{I;W^'''{^)), 
and (f)^''- converges to (p pointwise uniformly in each compact subset of $7 x [0, T]. 
From ()6.2|) we see that 

/ (cj(u^'=) + ek)V(t)'"^ ■Vipdx = for any if E H^{n). 
Jn 

By taking the limit A; — )■ oo, we obtain 

/ a{u)V(p • Vv7dx = 0, for any ip G H^i^) and a.e. t G (0,r). (6.17) 

Jn 

Therefore, for any function v G Cg°(r2), 

lim / V- ( (t>'"'{a{u'"') + ek)V(t)^>']vdx = - lim / (/)^'=(o-(u^'=) + efc)V0^'= • Vudx 
*:^°o Jn \ J >'^'^ Jn 

= — cl)a{u)V(j) ■ Vvdx 
Jn 

= - a{u)\7(p ■ [V{<t>v) - vV(p]dx 
Jn 

a{u)\V(j)\'^vdx. 



From dnU) we know that for any v G L°°{{0,T); C^(S7)), 



f r du^'' r r 

I / — — i'da;dt+ / / K{u^'=)e^yu'' -Vvdxdt 
Jo Jn ot Jo Jn 

/ V- {(j)^''—, ^V(l)''Avdxdt. 

By taking the limit fc — >• oo, we get 

/ / -^vdxdt+ I I K{u)Vu-Vvdxdt= / / cr(ti)|V0|\dxdt. (6.18) 

Jo Jn <^i Jo Jn Jo Jn 

Prom the regularity of u and (p, we know that the equations (j6.17p - (|6.18p actually hold 
for any <p G H^iQ) and v G L'^{iO,Ty,H^{n)). 

To conclude, we have proved the existence of a weak solution (u, cp) to the equations 
([ri])-([r3]) with the regularity I^M . 
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6.4 Uniqueness of solution 

Suppose that {ui,(j)i) and {u2,(p2) are two pairs of solutions to the initial-boundary value 
problem (jl.ip - (jl.3p . both satisfying (j3.6p . Let u = ui — U2 and (p = (/)i — cj)2- Then u and (j) 
are weak solutions to the equations 

du 

— - V • (K(ui)Vn) = V • {{k{ui) - k{u2))Vu2) 

+ {a{ui) - ct(u2))|V0i|2 + a{u2)V{(pi + (j)2) ■ V4> (6.19) 

- V • {a{ui)V4>) = V • ((a(ni) - a{u2))V(t)2) , (6.20) 
with the following boundary and initial conditions: 



u(x,t) = o, (/)(x,t) = o forxeao, tG[o,r], 

n(x, 0) = for X G il. 



(6.21) 



For any r G (0,T), we denote I-r = (0, r) and fi^ = fi x /,-. By applying Lemma 16.41 to 
the parabolic equation ()6.19p , we see that for q > n there exists 1 < p < oo such that 

\\u\\l^^q^) < C\\{k.{ui) - K,{u2))Vu2\\LP{U;L'}{n)) 

+ C\\i(^iui) - (T{u2))\V(l)l\'^\\LP(l^.Li/Hn)) 

+ C||a(u2)V((/)i + (1)2) ■ V^||ip(j^.i,/2(f^)) 

+ CT^^^P\\Vi(t>i + <P2)\\L^iU,L^in))\\'^4>\\L^PiIr;LHn)) 



where the constant C is independent of r. With the Holder regularity of ui, by applying the 
W^''^ estimates to ()6.20p . we obtain 



\L°°{ir;L^j{n)) < C\\{a{ui) -o-(-U2))V02||l°o(/^;L9(q)) < C||n||ioo(Q^). 
There exists Tq such that for r < Tq, the last two inequalities imply that 

ll'"llL°°(f7r) + \\^4'\\L°°{Ir;Li{n)) = 0. 

By dividing the interval (0, T) into small parts (r^jT^+i], A; = 0, 1, ... , each part satisfying 
Tk+i — Tk < Tq, we find that u{-,Tk) = ^(.,Tfc) = implies that u{-,t) = (j){-,t) = for 
t G [T/c,Tfc+i]. This proves the uniqueness of solution. 



7 Conclusions 

In this paper, we proved global existence and uniqueness of a weak solution to the degenerate 
thermistor problem by establishing a uniform-in-time BMO estimate for parabolic equations 
with possibly discontinuous coefficients. The physical hypothesis (H1)-(H2) are satisfied by 
metals and some semiconductors. The BMO estimate of parabolic equations established in 
this paper may be applied to many other equations of mathematical physics. 
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